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Abstract. In this paper we show that we can use a modified version of the h-p spec- 
tral element method proposed in |6 7 13 14] to solve elliptic problems with general 
boundary conditions to exponential accuracy on polygonal domains using nonconform- 
ing spectral element functions. A geometrical mesh is used in a neighbourhood of the 
corners. With this mesh we seek a solution which minimizes the sum of a weighted 
squared norm of the residuals in the partial differential equation and the squared norm 
of the residuals in the boundary conditions in fractional Sobolev spaces and enforce 
continuity by adding a term which measures the jump in the function and its derivatives 
at inter-element boundaries, in fractional Sobolev norms, to the functional being min- 
imized. In the neighbourhood of the corners, modified polar coordinates are used and 
a global coordinate system elsewhere. A stability estimate is derived for the functional 
which is minimized based on the regularity estimate in |2|. We examine how to par- 
allelize the method and show that the set of common boundary values consists of the 
values of the function at the corners of the polygonal domain. The method is faster than 
that proposed in 1 6 7 14 1 and the h-p finite element method and stronger error estimates 
are obtained. 

Keywords. Geometrical mesh; stability estimate; least-squares solution; 
precondi- doners; condition numbers; exponential accuracy. 



1. Introduction 

In B6I71 131 1411 h-p spectral element methods for solving elliptic boundary value problems 
on polygonal domains using parallel computers were proposed. For problems with Dirich- 
let boundary conditions the spectral element functions were nonconforming. For prob- 
lems with Neumann and mixed boundary conditions the spectral element functions had 
to be continuous at the vertices of the elements only. In this paper we propose a modified 
version of this method using nonconforming spectral element functions which works for 
general boundary conditions. 

For simplicity of exposition we restrict ourselves to scalar problems although the 
method applies to elliptic systems too. 

A method for obtaining a numerical solution to exponential accuracy for elliptic prob- 
lems with analytic coefficients posed on a curvilinear polygon whose boundary is piece- 
wise analytic with mixed Neumann and Dirichlet boundary conditions was first proposed 
by Babuska and Guo [3| within the framework of the finite element method. They were 
able to resolve the singularities which arise at the corners by using a geometrical mesh. 
This problem has also been examined by Karniadakis and Spencer in ifTTl . 
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We also use a geometrical mesh to solve the same class of problems to exponential 
accuracy using h-p spectral element methods. In a neighbourhood of the corners modified 
polar coordinates (t^, 9^) are used, where Zk = m r k an d ( r k, &k) are polar coordinates with 
the origin at the vertex A^. Away from sectoral neighbourhoods of the corners a global 
coordinate system is used consisting of (x\,X2) coordinates. 

We now seek a solution which minimizes the sum of the squares of a weighted squared 
norm of the residuals in the partial differential equation and the sum of the squares of the 
residuals in the boundary conditions in fractional Sobolev norms and enforce continuity 
by adding a term which measures the sum of the squares of the jump in the function 
and its derivatives in fractional Sobolev norms to the functional being minimized. These 
computations are done using modified polar coordinates in sectoral neighbourhoods of 
the corners and a global coordinate system elsewhere in the domain. The spectral element 
functions are nonconforming. For the modified version of the h-p spectral element method 
examined here a stability estimate is proved which is based on the regularity estimate of 
Babuska and Guo in [2|. The proof is much simpler than that of the stability estimate in 
M6I7I . Moreover the error estimates are stronger. 

The set of common boundary values for the numerical scheme consists of the values 
of the function at the vertices of the polygonal domain. Since the cardinality of the set of 
common boundary values is so small we can compute a nearly exact approximation to the 
Schur complement. Let M denote the number of corner layers and W denote the number of 
degrees of freedom in each independent variable of the spectral element functions, which 
are a tensor product of polynomials, and let W be proportional to M. Then the method is 
faster than the h-p spectral element method in M6I7I141 by a factor of 0(W 1 / 2 ) and faster 
than the h-p finite element method by a factor of 0(W). 

We now outline the contents of this paper. In §2 function spaces are defined and differ- 
entiability estimates are obtained. In §3 we state and prove stability estimates. In §4 the 
numerical scheme, which is based on these estimates, is described and in §5 error esti- 
mates are obtained. In §6 we examine the issues of parallelization and preconditioning. 
Finally §7 contains technical results which are needed to prove the stability theorem. 

2. Function spaces and differentiability estimates 

Let H be a curvilinear polygon with vertices A\,Ai,...,A p and corresponding sides 
Ti , T2 , . . . , r p where T, joins the points A,_ i and A, . We shall assume that the sides T, are 
analytic arcs, i.e. 

r ! - = {(<p^) ) ^))|§e7=[-i,i]} 

with <p,-(|) and \jfi(^) being analytic functions on 7 and |<p/(<i;)| 2 + | V / /(^)| 2 > a > 0. By 
r, we mean the open arc, i.e. the image of / = (—1,1). 

Let the angle subtended at A ; - be C0j. We shall denote the boundary dQ. of Q. by 

T. Further, let T = r[ ]|jr [1] , = \J ieS T u = Uie^i where @ is a sub- 
set of the set = 1, . . . ,p} and ,JY = = 1, . . . ,p} \ 2). Let x denote the vector 
x = (x u x 2 ). 

Let S£ be a strongly elliptic operator 

2 2 

-^( M ) = Y, ( a rA x ) u x,)x r + Y b r( x ) u x r + c{x)u, (2. 1) 

r,5-=l r=\ 
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where a hr (x) = a r , s (x),b r (x),c(x) are analytic functions on £2 and for any §1,^2 G K and 
any x E £2, 

£ /*>(£i 2 + & 2 ) (2-2) 

r..v=l 

with /Xq > 0. In this paper we shall consider the boundary value problem 
££u = f on £2, 
u = g® onrH 



dA/7 



onD'l (2.3) 



4 



where (du/dN)A denotes the usual conormal derivative which we shall now define. Let 
A denote the 2 x 2 matrix whose entries are given by 

A r<s {x) =a r , s {x) 

for r, s — 1,2. Let N = (A^A^) denote the outward normal to the curve r, for i G jY . 
Then ( is defined as follows: 



— x = L AU^. (2.4) 



Moreover let the bilinear form induced by the operator Jz? satisfy the inf-sup conditions. 
It shall be assumed that the given data / is analytic on £2 and g^ l \l = 0, 1 is analytic on 
every closed arc T; and g' ' is continuous on . 

By H m (£l) we denote the Sobolev space of functions with square integrable derivatives 
of order < m on £2 furnished with the norm 



|«» 2 - ^ "" a "" 2 



\a\<m 



Define r ; (x) to be the Euclidean distance between x and the vertex A, of £2. Let j3 = 
(/3i, j32, . . . ,P P ) denote a p-tuple of real numbers, < /3, < 1, i = 1, . . . ,p. For any integer 
k, let j3 + k = (J3i + it, j3 2 + k, . . . , fi p + k) . Further, we denote 

i=l i=l 

Let H™''(Q.),m > I > 0, / an integer, denote the completion of the set of all infinitely 
differentiable functions under the norm 

m 

K p v ; |a|=Jt,*=/ ' 

m 

V (Q) |a|=W=0 () 
For m = I = we shall write Hjj'°(£2) = Lp (£2). 



112 



P K Dutt, N Kishore Kumar and C S Upadhyay 



in- 1,1- 



I ,_I 



Let 7 be part of the boundary T of £2. Define H (y), m > 1,1 > to be the set of 



in. I 



all functions on y such that there exists / E H (Q.) with = f\y and 



For / an integer < I < 2, let 

V^(fl) = {ii(jf)|MGflp , '(fl),m > 

and 

s^(£2) = {u( x )\u 6 v^(fl), II Id^I^+wIIl^) < ca*-'(t-z)! 

for|a| = k — 1,1 + 1, . . ,;d > 1, C independent of k}. 

Let Q C R 2 be an open set with a piecewise analytic boundary <?2 an d 7 be part or 
whole of the boundary dQ. Finally 23 ^ 2 (y),0 < / < 2, denotes the space of all functions 
<p for which there exists / £ 23^ (2) such that / = (p on y. 

Next as in we introduce the space : 

<^(£2) = {u G //^ 2 (£2)| |fl a «(*)| < cAlf^^W)-', 
|a| =&= 1,2,...;C > l;d > 1 independentof £}. 

The relationship between and 23^ is given by Theorem 2.2 of |3j which can be stated 
as follows: 

<B2(£2)ce;2(£2). 

We need to state our regularity estimates in terms of local variables which are defined on 
a geometrical mesh imposed on £2 as in §5 of [3 1. £2 is first divided into subdomains. Thus 
we divide £2 into p subdomains S , ... ,S P , where S l denotes a domain which contains the 
vertex A, and no other, and on each S' we define a geometrical mesh. Let 6* = {£2* -,_/' = 

1, . . . i = 1 , . . . ,4.;} be a partition of S k and let 6 = Uf=i ® ■ The geometrical mesh 
imposed on £2 is as shown in figure 1 . 

We now put some restrictions on 6. Let (r^, 9%) denote polar coordinates with center 
at Afc. Let = lnr^. Choose p so that the curvilinear sector £2^ with sides and T^+i 
bounded by the circular arc B k p , center at Aj and radius p satisfies 



a k c (J £2* ; , 

£2 A may be represented as 



Q. k = {(xi,x 2 )e£i:0<r k <p }. 



(2.5) 



Nonconforming h-p spectral element methods 



113 



P r. SI 3 A > 




Figure 1. Geometrical mesh with M layers in the radial direction. 



Let tfj ;{ , 1 < / < 4 be the side of the quadrilateral Qjj G 6. Then it is assumed that 



1<<?<1, 1 = 1,3 



(2.6a) 



(x l =hl J (pl Jl (r]), 

tx- " -i<n<h 1 = 2,4 

and that for some C > 1 and L > 1 independent of i,j,k and /, 

<CZ/f!, f = 1,2, 



^ it ^ 




d ' k ,s 




i 





(2.6b) 



(2.7) 



We shall place further restrictions on the geometric mesh imposed on Q. k later. Some of 
the elements may be curvilinear triangles. 

Let (r k , 9 k ) be polar coordinates with center at A k . Then £l k is the open set bounded 
by the curvilinear arcs F k , T k+ i and a portion of the circle r k = p. We divide Q. k into 






e k =^(r k ) 























Figure 2. Curvilinear sectors. 
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curvilinear rectangles by drawing M circular arcs r k = aj' = p pL^ + J ,j = 2, . . . ,M+ 1, 
where ji k < 1 and I k — 1 analytic curves C2, ■ ■ ■ , Cj k whose exact form shall be prescribed 
in what follows. Let a* = 0. Thus I k / = I k for j < M; in fact, we shall let 4 j = 4 f° r 
j < M+ 1. Moreover 4j < / for all k, j where / is a fixed constant. Let 

T k+j = {(r k , e k )\d k = ff(r k ), 0<r k < p}, 

j = 0, 1 in a neighbourhood of Aj- in £l k . Then the mapping 

» = ft, h = - Vf)7f (ft) - (4 - VM(ft)], (2-8) 

where /* is analytic in p,t for j = 0, 1, maps locally the cone 

{(p*, fo): < Pi < <r, y/f < <fe < v/*} 

onto a set containing £l k as in §3 of 0. The functions fj satisfy/^ (0) = y/f, ,/f (0) = 

and (/f)'(O) = for 7 = 0, 1. It is easy to see that the mapping defined in ( 12. 8t has two 
bounded derivatives in a neighbourhood of the origin which contains the closure of the 
open set 

& = {{pk, <h)- 0<Pk<P>Vi<h< wt)- 
We choose the 4 — 1 curves C2,--,Ci k as 

Q:<j) k {r k ,d k ) = \l4 
for i = 2, . . . ,4- Here 

Let Ay/f = y^, j - y/f. Then {vf}i,A are chosen so that 

max(Ay/f) < A(min(Ay/f)) (2.9) 

i,k i,k 

for some constant X. Another set of local variables (t k , 9 k ) is needed in a neighbourhood 
of Q. k where 

T k = \nr k . 

In addition, we need one final set of local variables (v k ,(j) k ) in the cone 

{{pkAkY o<p k <p,wi<h< 

where 

Let = {(r k , 9 k ): < r k < ji} fl£2. Then the image 5^ in (v^, variables of is given 

by 

Si = {(VkAk): -°°<v k < Inji.yf < 0* < 
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Now the relationship between the variables (Tj, %) and {Vki^k) is given by (^,8^) = 
M k (v k ,<j> k ),viz. 

% = Vjt, 

o k = i-f^ir^ vf)/i> v *) ~ (<k - vft/SW]- (2-10) 

Hence it is easy to see that J (v k , § k ), the Jacobian of the above transformation, satisfies 
Ci < \J k (v k , (j> k )\ < C 2 for all (v t) fc) e 5* , for all < /x < p. 

We now need the fundamental regularity result from 0, viz. Theorem 2.1 which we 
state as follows: 

Iff £ fl^°(n), e H™ + ^ 3 '^ 3 {T^), j = 0, 1, < ft < 1, ft > ft* and m > 0, then 
the solution of (2.3) exists in Hp +2 ' 2 (Q.) and 

l 



[./: | 

f> l " J \ "13 v> j=0 h"" 2 J ' 2 J (rW) 



Let us define a, = 1 — ft . 

We now state the differentiability estimates for the solution u of (12. 3t which will be 
needed in this paper. 

PROPOSITION 2.1. 

Let 1 — a k > 0. Then for X k < Cft, 



£ |D £ v ;^( M - M (A,))| 2 e- 2 ^^dv,d0, 



V/ - 7 - 00 |e|<m 

< li 2 »(CV"- 2 (m-2)!) 2 (2.11) 

for < jl < p with y k < OC k - X k .Ifl — OC k <0 then for X k < 1 /2, (2.11) remains valid for 
< /I < p w?/z % = 1 /2. 

The proposition can be proved in the same way as Theorem 2. 1 of |6|. 

3. The stability estimate 

Let 

2 2 

-^( M ) = - £ (ai,jUx } )xi + Y,biii Xi + cu (3.1) 

U=i '=1 

be a strongly elliptic operator. We now consider the following mixed boundary value 
problem: 

ifw = / in Q., 

%u = u\ r [ ] = and 
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Here the conormal derivative y\ u is defined as follows. Let Ji C rW and let N = {N\ , A^) r 
denote the unit outward normal at a point on Ji- Then 

fdu\ 2 

7i" = (^J a = £ (3.3) 

Moreover, let the bilinear form induced by the operator _£f satisfy the inf-sup conditions. 
We can now state the regularity result Theorem 2.1 of J2 as follows: 
Let u be the solution to (3.2). Then 

The above estimate for £ = is used to prove the stability estimate Theorem 3.1. 

We remark that in Theorem 5.2 of (9), Guo and Babuska have extended the above 
regularity result to elliptic systems. Hence the method applies to elliptic systems too. 

Divide the polygonal domain £2 into p sectors Q. 1 ,£l 2 , . . . ,£1'' and a remaining portion 
Q. p+l . Further divide each of these subdomains into still smaller elements 

{Oy, l<i<hj, i<j<M,i<k< P }. 

Let 

£l p+l = {Cllf. \<k<p,M<j<J kl \<i< I kJ }. 

We shall relabel the elements of Q. p+1 and write 

Q. p+l = {Q.f +l :l<l<L}. 

Now define the space of spectral element functions Yl M - w = /(v^, <j>ic)}ij,k.i 
{u^ +1 (4 ,17)}/}, where u k ix = h k a constant for all i and 

Wj Wj 

Uij(Vk, <Pk) = £ L *v v k #> 1 < i < M - 

r=li=l 

Here 1 < Wj < W. Moreover there is an analytic mapping Mf +l from the master square 
S = (-1, l) 2 to £2f +1 . We define 

w w 

«f +1 K +1 (^,T])) = ££^^^. 

r=ls=l 

Let w G ff*' 2 (£2). Now for 1 < j < M, 

k r 2 / k \J?w\ 2 dx = / r 2( - 1+ft ' ) |JfV| 2 dT A .d0 i . (3.5) 

Here £2f ■ is the image of Q. 1 - • in (tj, 0^) coordinates and _Sf*w = r^Jzfw. It has been shown 
in that if we let y 1 = and 3^2 = 0^ then 

^=-| i |:(4^) + ES!» J ,+^. (36) 
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Let O denote the matrix 



O k 



cos 6^ — sin Ok 
sin Ok cos Ok 



and 



A k = 



^2,1 ^2,2. 



Then A* = (£>*) r A0* 



Let J k {v k ,$ k ) denote the Jacobian of the map M k (Vk, <j> k ) defined in §2. Then for 1 < 



(3.7) 



Here Clfj is the image of Q. k j in (v*, 0t) variables and 



Now 

J? k jw(v k ,$ k ) = A k jWVkVk + 2B k jWvkh + C k jW Mk 
+ D k j w Vk +E k J w h +F k j w. 

Let A k j be the polynomial approximation of A* of degree W ; - in and 0^ separately, 
as defined in Theorem 4.46 of [ 12 1. Now we define a differential operator with polynomial 
coefficients (-2f ,-) , which is an approximation to J?f ; - as follows: 

{^ k jY'w = A k jWvkVk +2BijW Vkt f, k +C k jW Mk +D k jW Vk 
+ E k jW^ k +F k jW . 
Let X k = 1 - jSjfc. Then for 1 < j <M, 



|^Wv*,0*)|V 2 ** v *dv*d0* 



- f \(^j) a w(v k ,<k)\ 2 ^ UkVk dv k d^ 



<e w ((PMr +1 ^r 24 l|w(v,,#)-w(A,)||2. t 
+ (p M f +1 -0 4 - 2 ^|w(A,)| 2 ). 



(3.8a) 



Here e w — > as and, in fact, e w is exponentially small in W. 
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Moreover, if w(Vk, tyk) =w(A^), a constant in £2*j for 1 < z < 4> then 
£ / |^Xv,,^)| 2 e- 2 ^dv,d^<e M |w(A,)| 2 . 



Here e M — > as M — > °° and e M is exponentially small in M. 

Hence we conclude that if w(v^, ^) = w(A jt ), a constant in £l k n for 1 < i < 4, then 

EE L |^>fj(v fc ,#)| 2 e- 2 ^dv fc d^ 

/4m \ 

[LL(p^ +1 ~ J r 2 ^tj-w(A k )f +HA k )\ 2 ) 

\i=lj=2 2 - n -.J ) 

+ e M |w(A,)| 2 . (3.8b) 



Here C is a constant. 
Now 



/ , |j^H 2 dxidx 2 = f |^wf +1 | 2 7f +1 d^d77. 

Jof +1 Js 



Here Jf + 1 (| , 77 ) is the Jacobian of the mapping Aff + 1 from 5 to Of + 1 . Let Jzlf + 1 (| , 77 ) = 

Jz?(<i; ,17) \J Jf +l . Once more we can define (Jz?/ ,+1 ) a , a differential operator which is an 

approximation to Jz?/^ 1 in which the coefficients of -5f/' +1 are replaced by polynomial 
approximations. It can be shown as before that 



E L, I^H 2 d* ld * 2 

;=1 ,/ "; 

ct / i(^/ ,+i ) n < +i (^T7)i 2 d^dT 7+ew x:ii< +i (^T 7 



CT / ( '/'!' ')V '(£ ( ]) 2 dcdn + f„ ri»rr '(c.r))|| 2 s . 

Here C is a constant and £„, — > as W — In fact, e w is exponentially small in W. 
We now prove a result which we shall need in the sequel. 

Lemma 3.1. Let (0 € Hp' 2 (Q.). Then there exists a constant C such that 
iff (\co(A k )\ 2 + £ / |D« A (a)(v,,«-»(A fc ))| 2 

xe-^dv^j+HcoCx^xz)!! 2 ) 
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(3.9) 



Here X k — 1 — fa. 

Proof. Let y k e Cq(X) such that \)f k (r k ) = 1 for < p and l^fo) = for r k > p 1 for 
fe= 1,2,... , p. Here p 1 > p is chosen so that t = {(jci,Jt2): r i _ p 1 } have the property 
that £l k pl nfl', = if & 7^ /. We define co k = CO ify for k = 1, . . . ,p and C0q = 1 - ££ =1 Cty. 
Then to* € fl^(£i) for k = 1, . . . ,p. 

. . 2 2 

Now by Lemma 2.1 of |3J, (£2) C C(£2) with continuous injection. Hence we con- 
clude that 



ii^(A,)i 2 <ci:iito A ii 2 




k=l k=\ 



Therefore 




(3.10) 



We now cite Lemma 2.2 of 10). Let u £ H' (£1). Then 



(i) 



£ WD^u^W <C\\u\\ 



or 1=1 



(ii) Let u(Ai) — 0, for z = 1 , . . . , p. Then 



l«^- 2 IU <c||«|| 



Hp (a) 



From (i) we obtain 



i L e i^ A 0(v,,^)iv 2 ( i -^dv^<c(iifi)ii 2 22 ). (3.H) 

*=l- /n ^ |oe|=l V (£i) 



Here C is a generic constant. Now using (ii) we get 



/ \co k {v k ^ k )-co{A k ) Wk \ 2 ^ l -^Av k d^ k 
Jo. 

<C(||0),|| 2 2 2 +\co(A k )\ 2 ). 
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Hence 



Finally, 



p 

E 

k=\ 



P 

E 

k=l 



a* 



\co(v k ,<j) k )-a){A k )\ z e 



2 p "2(l-ft)v t 



dv^d^ < C\\co\ 



«p 2 (&) 



£ |Z)« A co(v,,fe)| 2 e 



2_-2(l-ft)v t 



dv t dfc<C||a>| 



a\=2 



2 2 
Hi' («) 



(3.12) 



(3.13) 



Combining the estimates (3. 10)— (3. 13) we get (3.9). 



We now introduce some notation which is needed to state the stability estimate Theo- 
rem 3.1 which is the main result of this section. 

Let Y s be a side common to the elements £2m +1 and Q.„ +1 and let y s C Q. p+1 . We may 

p+i 



assume that % is the image of tj = — 1 under the mapping M„ + 1 which maps 5 to £2m +1 and 
also the image of 77 = 1 under the mapping M„ + l which maps S to Q.„ +l . By the chain 
rule 

K +l ) xl =K +1 k^+« +l )rin Xl , and 

Now let gx } denote the polynomial approximation of (§,Tj), of degree W in % and 77 
separately, as defined in Theorem 4.46 of fPH . In the same way f) Xl ,£, X2 and fj v , can be 
defined. We now define 

= K + \ L + « + % f\M , and 



Let 



= \K +i ^,-i)- u p +l ^,i)\\ 



0.;' 



,P+l\a]\\2 



[LI 



J "1/2, R 
x 2 ' Jll l/2, rs 



, - 1) - (i^- 1 )- , 1) llf 



I 2 

1 1/2,/ 



and 



)^,-i)-K +1 )^,i)||^ ; 



I 2 

n/2,1 



Here / = (-1, 1). Next, let 7, C n <?i2 p+1 and let % be the image of 7] = -1 under 
the mapping Mf, +1 which maps 5 to iim +1 . We can define ( ^g^ )°, an approximation to 



~3T 



as before. Let 



ll« P+1 || 2 + 



duP +l \ a 
dT 



In the same way, if y s C n <3£F +1 , 



1A» 

/ dum +1 \ 

[dT ) 

I ( dl ' p+> V I 

I \ dN )a\ 



-1) 



l/2>» 



1/2,/ 

can be defined. 
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Let Ys C r' 1 ' (~l dQ k for 1 < k < p. Let % be the image of 7 in (t^, 9k) coordinates and 
% be the image of y s in (Vk, fa) coordinates. Let (n\ , «2) be the normal at a point p on %. 
Define 



(du k \ X k du k 



Now % is a portion of the straight line fa = a, where a is a constant. Let (^)^n, denote 



an approximation to (-fc)^ as before, and using this || || . can be defined. Let 

y s C H*\ Define 



Let 



d(A t ,%) = inf {distance^,*)}. 



P M 4 



E E E(p< +1 ^)^iK^) fl «f,(v„^)ii 2 

k=\j=2i=\ 'J 



E E 



-2A, 



/t=l 



r ! ca*uB p ,M(r s )< c > 



.fcll|2 



x(iiKjii^ + ii[«r]ii^ 



+ E E (N 2 + E ^a*,/,)- 2 ** 



x (||^-^|| 2 , + ||<|| 2 Afe )) 

+ E E E 

/e^ *=/-! rt c5n*nr,,^(y t )<oo 



(3.14) 



Here {{«^.(v t , <MH M , K + ' (§, J?)}/} G n M ^ w and ^ = % for 1 < i < I k . Moreover 
jtt(7i) denotes the measure of %. Next, we define 



^ Mi^,fa)}i d M + \^^)}i) 

=i\\(^r i ) a uf +i ^,r ] )\\i 



1 1 0.5' 



+ E (H[»^ 1 ]t + ii[« 1 ) fl ]iiU. + ii[« 1 ) fl ]ii? A ,) 
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I I 

E E 



u 



P+l||2 



Let 



"0,J» 

du" +1 Y 

A 

P+i , 



dT 



1/2,14, 



dN 



1/2* 



r ' (K;(v,,^)}u, ft ,K +1 (§,r?)}/) 
= 'CI({»t-(v fe ,^)} ! - M ,K +1 (^T,)} ; ) 

+ CI(R'( V ^^)}^,K +1 (^^)}/)- 

We can now state the main result of this section. 
Theorem 3.2. For M and W large enough the estimate 



h M 

E [ i^i 2 +EE(Pi"f +1 ^r 24 ii«f.. ; (v,,^)-/^n 



I 

k=\ 



i=l;=2 



2,n* . 



L 

I 

/=! 



< c(in w) 2 r* w ({i4j(v k , & >*»)}/) 

holds. Here C is a constant. 



(3.15) 



(3.16) 



Proof. By Lemma 7.1 there exist {{v; ; (Vi, <j>k)}ij,ki { v 1 (^j 1 ?)}/} sucn that w defined 
as w — u + v € /^' 2 (£2). Moreover v*j = for all 2 and A;. Hence by Theorem 2. 1 of 0, 



lk|| 2 22 <C 



V 



Now vf j (v*, 0*) = for 1 < i < l k . Hence by (3.8), 



dw 
~dN 



(3.17) 



||J2? W || 



/ p M h 

<2 III ii^tjHMM' +E (^)ll 

U=U=2/=1 



J 1=1 



2 

o,s 
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(P 4 M L 
III(pMf +1 - 7 r n iv*,(v fc ,^)ip +ek +1 (^)ii 2 2 . 
k=\i=\)=2 'hi 1=1 

(P M 4 P \ 

E E E(pMf +1 -r n *K/v*,<« -m 2 + E N 2 
Ar=l y=2i=l hi k=l ) 



Now using Lemma 7.2, 



|W|| 2 3 ; 



ff 2 ' 3 (r[°l) 



dw 
dN 



fl|'2(rW) 



<C(lnW) 2 ( £ N 2 +£ £ 

r s £3n*nr,,^(y s )<oo 



+ E E E d(A k , 7s )- nk 

leJ^k=l-\ 7s cdQ.i<nr,,Li{y s )<<» 



du' 
dn 



k\ a 



A* 



1/2,* 



+ E E 



' e ® 7s can p+1 nr ; 



i/2,rs/ 



E E 

,e ^ ri c3n p+1 nr, 



A 



dN 



x(ii[/]ii^ + ii[«r]ii^ fe + 



l/2,» 

24 



y s C£2' 



E (Il^lt + IIKOIII^ + IIKOIIIU,)) 



(p p h M 

Ei^i 2 +EEE(pMf +1 " y r 2 ^ll^(v„« 
fc=l k=li=lj=2 



+in«r i «.i?)ie, ■ 



124 P K Dutt, N Kishore Kumar and C S Upadhyay 

Combining (3. 17)— (3. 19) we obtain 



( P h M L \ 

+c III(PMf +1 - y ')- 2A ivt,(v fc ,^)nj +EK +1 (^)iiL 

\k=li=lj=2 hi 1=1 J 

( p M ' k „ , ■ p \ 

+M L L E (ppf + J y 2h \\ u\j {v k ,(j) k )-h k \\ 2 + e N 2 

\<r=l ;'=2i=l 'hi k=l ) 



Now using (3.9), 



p / 4 m \ 

E i^i 2 +EE(p< + l - J r 2 H4j(vM-h\\ 2 

k=\ \ i=l ;=2 2 '"'j/ 

+EK +1 (^)ll 2 , 

/ P h M 

< * HI* + E E E(p< +1 ")- 2A ivf j( v^,)iij^ 

\ <") *=li=l j=2 hi 

+ thf + \^r 1 )\\ 2 2 )- (3.21) 



Combining (3.20) and (3.21) gives 



p ( 4 m \ L 

E l^l 2 +EE(PMf +w )- 2 ^ll«?,(v,,^)-^|p +EK +1 (^)ii 2 s . 

k=\ \ i=lj=2 hi I 1=1 



< 



C(lnW) 2 V M ' W ({ul/VkAk^K 1 (l.n)}/) 

(P h M L \ 

E E E(ppf +w )- 2A i^(v^,)n 2 +e K +1 (^)iiL 
k=li=lj=2 'Ui l=l J 

(P M h P \ 

E E Lip^T'H^M-hkii 2 + E W 2 
A'=l ./=2;=1 hi k=l ) 

+% (ejk +1 (^)II 2 J + £ m (J>n • (3 - 22) 
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Now by Lemma 7.1, 

p 4 M L 

III(PMf +1 - y )- 2A ivf j( v^ fc )ii 2 +EK +1 (^)ii 2 s 

k=\i=lj=2 • U 1=1 

\*=ir ! cn*uB^. 1 u(f ! )<~ 

x(ii[/]ii 2 , fe + ii[(<nii 2 A , + iiHr]ii^) 



+ E (iiK +1 ]iiL + ii[(< +1 )iii 2 /2 , K! + ii[K +1 r]ii 2 /2 , X! ) 

/ P 4 M 



+% ( E E E(P< +1 ~T 2 *i«L-(v,,<« -m 2 

\Ar=ll=l;=2 hj 

+ thi +l tt^)\\\)j- (3-23) 
Combining (3.22) and (3.23) we get the result. 



4. The numerical scheme 

As in §3, 

= {(v k ,h)- V/ <v k < v* +1 , y/f < # < y/f +1 } 

for 1 < j < M, 1 < ; < 4j, 1 < < p in v,t and fa variables. 

We now define a nonconforming spectral element representation on each of these sub- 
domains as follows: 

u lj(Vk, fa) =h, if 7 = 1, 1 < i <h, 1 < k < p 

and 



uij(v k ,fa) = £ E a ^» v rc 



OT=1 W=l 

for 1 < j <M, \ <i < I kJ , \<k<p. Here 1 < W, < W. Let 

= {af + \i<i<L}. 

We define the analytic map Mf + from the master square S = ( — 1 , 1 ) 2 to 1 and let 

w IV 

< + Vf +1 (<?, 77)) = EE «»v>rv. 

m=l n=l 
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Let/f +1 (£,T,) 



f{Xf + \£,,T)),Yf + \$,T])) for 1 < I < L and Jf +l £,r}) denote the 



Jacobian of the mapping Mf + . Define F/'^ 1 ,77) = f l 1 >< ,, ; 

let F ; p+1 (<^,77) denote the unique polynomial which is the orthogonal projection of 
Ff^i^n) int o the space of polynomials of degree 2W in E, and 77 with respect to the 
usual inner product in H 2 (S). 

Next, let the vertex A k — (x k ,yk)- As defined in §2 we have the following relationship 
between (t^, 9k): and (\>k, <j>k) coordinates: 



J', ,+1 ^,ri) and 



Vk = X*, 



9k 



iwi-wt) 



l(h-v4)fi(e Vk )-(h-v!:,)ri(e 



Vk) 



Define f k (x h 6 k ) = e 2z *f(x k + ^ cos k ,y k + e T * sin 9 k ) for I <k<p, and F k j{v k , <M = 
f k (%k, 6k) for (Vk, <pk) G ^fj- Let F k j(Vk,(l>k) denote the polynomial of degree 2Wj in v,t 
and <j) k variables which is the orthogonal projection of F k j(v k , <j) k ) into the space of poly- 
nomials of degree 2Wj in \>k and <j) k variables with respect to the usual inner product in 
H 2 (Q. k j). Hem 2 < j <M. 

We now consider the boundary condition u = gk on F k forked and let = gk 

on r,t for k G JY ' . Define 



u = g k (x k + e v *cos(f k (e v *)),y k + e v *sm(f k (e v ><))), forked, 

— \ — „v ko ,J r ,_ l _ P v krnQ (fkf v k \\ „. j_„v t ^„ffki 
dnJ Ak 

for A: e 



= e v *gt(^ + e v *cos(/*(e v *)) , M + e v * sin(/*(e v *))), 



Let Z [ ; -(Vfc) be the orthogonal projection of l\ (Vk) into the space of polynomials of degree 
2Wj with respect to the usual inner product on H 2 (v k , v| +1 ) for 2 < j < M. 
Consider the boundary condition u = gk on T k n d£l k ~ l . Define 



i k 2(y k - l ) = { 



ft (x,_ 1 +e^-icos(/f- 1 (e v *-')),^ 1 +e v *-'sin(/f- 1 (e v ^))) 



for k e 9 

du 
dn 



■k-l( p v k -i\ 



= e^g^ +e^-i costf' 1 ^)) , >', 



+e v *-i sin^f- 1 (e v *-i ))) for e 



Let a,t = u(Ak) if 7*: or 7^ +1 e f^. We define /^(v^-i) to be the orthogonal projection 
of l k (Vk-i) into the space of polynomials of degree 2Wj with respect to the usual inner 
product on H 2 (v k " 



• k ~ l ,V k ~\)for2<j <M. 



Finally, let r^H^f = Q be the image of the mapping Mf 1 " of 5 onto Of 
corresponding to the side E, = -1. Let of(rj) = g*(X/' +1 (-l,T7 i ),y/' +1 (-l,?7)), where 
— 1 < TJ < 1. Define 0^(77) to be the polynomial of degree 2W which is the orthogonal 
projection of 0^(77) with respect to the usual inner product in H 2 {— 1, 1). 
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Now we formulate the numerical scheme for problems with mixed boundary condi- 
tions. 

Let {{vy(Vfc, <j>k)}i,j,ki {v P [ + ! S n M ' W , the space of spectral element functions. 
Define the functional 

P M h 

= E E I(PMf +1 - J r 2h K^jHj(y k ,h) 

k=l j=2i=\ " M iJ 

+f e </(A^r 2 ^ 
x(ii[v»]iii t +ii[(<)iiif At +ii[(v* fc niif At ) 

m 

+ E E I rf(A fc , 7v )- 24 (ii^-^) 

mG®<c=m-i rt c3a*nr m ,^(j> t )«» 

- (lZ-k+l- a k)\\l % + \\v k Vk - (l'm- k + 1 )vj 1 2 /2 ^) 
m m 

+ E E (^-"*) 2 + E E 

mE@k=m— 1 mG- 4' k=m—l 



r,c,9a*nr m ,fi(f J )<~ 



In the above jlt(y s ) denotes the measure of %. 
Next, define 
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Let 

+ CljHj^ («,»»)}»)• (4-3) 

We choose as our approximate solution the unique {{zf ;(V£,0£)},-j i fc,{z[ +I (4, 17)}/} S 
n M W ', the space of spectral element functions, which minimizes the functional 

A brief description of the solution procedure is now given; a more detailed examination 
is provided in §6. The above method is essentially a least-squares method and the solu- 
tion can be obtained by using preconditioned conjugate gradient techniques (PCGM) to 
solve the normal equations. To be able to do so we must be able to compute the residuals 
in the normal equations inexpensively. In B5I141 it has been shown how to compute these 
efficiently on a distributed memory parallel computer, without having to filter the coef- 
ficients of the differential operator and the data. The evaluation of the residuals on each 
element requires the interchange of boundary values between neighbouring elements. 

The values of the spectral element functions at the vertices of the polygonal domain 
constitute the set of common boundary values Ub- Since the dimension of the set of com- 
mon boundary values is so small a nearly exact approximation to the Schur Complement 
matrix can be computed. Now on the subspace of spectral element functions which van- 
ish at the set of common boundary values it is possible to define a preconditioner for the 
matrix in the normal equations such that the condition number of the preconditioned sys- 
tem is (9 ((In W) 2 ). Moreover, the preconditioner is a block diagonal matrix such that each 
diagonal block corresponds to a different element, and so can be easily inverted. 

Hence an exponentially accurate approximation §" to the Schur Complement matrix § 
can be computed using 0(W In W) iterations of the PCGM. To solve the normal equations 
the residual in the equations for the Schur Complement SUb = hg must be computed to 
exponential accuracy and this can be done using (9(WlnW) iterations of the PCGM. The 
common boundary values Ub are then given by Ub = (E> a )~ 1 hB- The remaining values can 
then be obtained using <9(WlnW) iterations of the PCGM. 

5. Error estimates 

Let {{^j(v k ,fa)}i Jtk ,{^ +l (^,r])}i} minimize tr"'* ({v£ ; .(v t , &)},-, jJc , {vf 1 (|, 77)},) 
over all {{vf j(v k , fa)}ijjc, { v f +1 (4 1 T)}/} £ n M W , the space of spectral element func- 
tions. Here z k n = b k for all i, i(v k ,fa) is a polynomial in v k and fa of degree Wj, 

Wj < W and (§,tj) is a polynomial in | and rj of degree W as defined in §3. We 
choose W proportional to M. Then we have the following error estimate. 

Theorem5.1. Let a k = u(A k ). Let Ufj{y k ,fa) = u(v k ,fa) for (v k ,fa) G A*- and 

U'i +l (E,,r\) ) for (§,Tj) eS. Let aj < Wj < W for some positive a for j > 2. 

Then there exists positive constants C and b such that for W large enough the estimate 

i\b k -a k \ 2 +iH(p^ +i -r 2h 

k=\ k=\ ;=2i=l 
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x\\(4j-U!j)(v k ,fa)-(b k -a k )\\l^ 

+ t\\tf +1 -Ur l )($M\ls<Ce- bw (5.1) 
1=1 

holds. 

We use the differentiability estimates stated in Proposition 2.1 to prove the result. The 
proof of the above Theorem is very similar to the proof of Theorem 3. 1 in [ 14 1 and hence 
is omitted. 

Remark. We can construct a set of corrections {{cf :{Vk,$k)}i.j,k, (I i S n M,H ', 
the set of spectral element functions, so that corrected solution {{^j(Vk,<j>k)}i,j,k> 
77)},} defined by 

={{^(^^)}w,{«r 1 («^)M+{{to»^)w{ c r 1 (^ T J)}/} 

is conforming and belongs to H 1 (CI). These corrections are defined in §3.5 of [ 14|. Then 
the error estimate 

\\(u-z)(^y)h.n<Ce- bw 

holds for W large enough. Here C and b denote constants. These constructions are similar 
to Lemma 4.57 in ifTIl . 



6. Parallelization and preconditioning 

Let U be a vector assembled from {gyfc}f =1 , where u k n = gk for all i, and the val- 
ues of {{u\ /(Vfc,^)}jj£,{Hj' +1 (£;,Tj)},} at the Gauss-Lobatto-Legendre points 
are arranged in lexicographic order for 1 < k < p, 2 < j < J^, 1 < i < 4;. Let 

{Uij(v k , h)}i,j,k, {z"i +[ (4 , rj)}/} minimize t M ' W ({v^(v fc , h)}i,j,k, K + ' (4 ■ »?)}/) over 
all {{v^j(Vk, (j>k)}ij,k, { v 'l +> (4 7 1 7)}/} G n M ' W , the space of spectral element functions. 

Let Ub denote the values {gfc}f =1 and Uj the remaining values of U. We now define a 
quadratic form 

^' M, ({«f. ; (n,^)} 1 ,M,{«r i (4,T7)}/) 

P P M 4 

= E bl 2 + E E E(PMf +1 - ; )- 2A i«f, / (4,r 7 ) -^ni iS 

k=\ k=lj=2i=l 

+ Ell»r +1 (4,r])|!l 5 . (6.1) 

/=i 

It should be noted that u k j (Vk, (j>k) — gk for 1 < i < Ik- Moreover for /' < M, % is a linear 
function of \>k and tj is a linear function of (j>k such that the linear mapping Mfj(^,rj) 
maps the master square S onto Q. k ■. 
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To solve the minimization problem we have to solve a system of equations of the form 
AZ = h. (6.2) 
Here A is a symmetric positive definite matrix and 



where f M,W {{u k i;j {v k , <j> k )}ij, k, K + ' (£,??)}/) is as defined in (3.15) in §3. 
Now A has the form 



(6.3) 



A = 



Au A IB 
Abi Abb 



(6.4) 



corresponding to the decomposition of U as 
U 



Ui 
U B 



and h has the form 

h = 



hi 



To solve the matrix equation (6.2) we use the block L-U factorization of A, viz. 
/ C 



A T A' 1 



An 

§ 



/ A n l A IB 
/ 



where the Schur Complement matrix § is defined as 

E>=Abb-AJbA h 1 A ib . 



(6.5) 



(6.6) 



To solve the matrix equation (6.2) based on the L-U factorization of A given in (6.5) 
reduces to solving the system of equations 



where 



§Z B = flB, 



hg = Iib —AjgAjjhj. 



(6.7) 



(6.8) 



The feasibility of such a process depends on our being able to compute AjbVb, AjjVj and 
AbbVb for any Vj, Vb efficiently and this can always be done since AV can be computed 
inexpensively as explained in ch. 3 of Tl4\ . 

However in addition to this it is imperative that we should be able to construct effective 
preconditioners for the matrix A// so that the condition number of the preconditioned 
system is as small as possible. If this can be done then it will be possible to compute A^'V/ 
efficiently using the preconditioned conjugate gradient method (PCGM) for any vector V/. 

Consider the space of spectral element functions Il ' , such that for {{w* (Vfc, <j>k)}i,j,kj 
{u^ +1 , ^7)}/} £ Tl^' W we have k*j = for all i and k. Let U be the vector corresponding 
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to the spectral element function {{wfj(V£, &)};,./,*> {"/ (4 , Then Ub — and t/ = 

and so 



y (K, 7 (v,,^)} (iM ,{«r (4,T?)} / ) = c//A // f/ / . (6.9) 

Now using Theorem 3. 1 we have the following result. 

Let {{«* J (v fc ,#)} ( - M ,{«f 1 (4,T/)},} enf W . Then the estimate 

P M h L 

III(PMf +w )- 2A l4/^n)lli, s +III"f +1 (4,r I )||!, s 

fc=ly=2i=l (=1 

< C(ln W) 2 r MM ' ({«fj(Yt, t )}i,M. K + 1 (4, T?)}/) (6- 10) 

holds for W large enough. In the above, u k n — for 1 < k < p and 1 < i < 4. 
Let us define the quadratic form 

^ M ' W ({4 y -(V fc ,#)} ( - M ,{«f 1 (4,T ? )} ; ) 

P M h L 

= LLL(p«!' +w r 2a *ii«fj«.»?)iii s +Lii«r 1 «,»?)ii| s (6.id 

k=\j=2i=\ 1=1 

forall{{M^.(v,,^)} I - M ,K +1 (4,Tj)} / }Gn^ w . 

Now using the trace theorems for Sobolev spaces it can be concluded that there exists 
a constant K such that 

^({ M f J .(v„^)} ( - iM ,K +1 (4,T / )} / ) 

<^^' W (K/v, ) ^)} i ', M ,K +1 (4,T?)}/) (6.12) 

for {{«* y (v 4 , {«r ' (4 , r?)}/} € I# w . 

Hence using (6.10) and (6.12) it follows that there exists a constant C such that 

ir M ' H/ ({«f J (v,,^)} 1J .,,K +1 (4,7 7 )}/) 

<^ M ' w, ({«f J (v,,^)} IJ> ,K +1 (4,T 7 )} / ) 

<C(lnW) 2 ^ M ' M '({«f J (n,^)} 1 ,M,K +1 (4,'7)}/) (6.13) 

for all {{«f J (v,,^)}, M ,K +1 (4,7 7 )}/} 

Thus the two forms f M ' W ({«* y (v t , {"f (4 , 1?)}/) and ^ ({^/v fc , A !| , ;i . 

{f^ +1 (4 are spectrally equivalent. 

We can now use the quadratic form fy"'™ ; -(v^, <j>k)}i,j.k, { M / + ' (4 1 *?)}/) which 
consists of a decoupled set of quadratic forms on each element as a preconditioner 
for A//. This can be done by inverting the block diagonal matrix representation for 

^({«f J (v,,^)} I - M ,K +1 (4,7 7 )} / ). 
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Now from (6.13) we can conclude that if we were to compute (Ai/)~ l Uj using the 
PCGM then the condition number of the preconditioned matrix would be 0((\nW) 2 ). 
Hence, to compute (A//) _1 f// to an accuracy of 0(e~ bw ) would require 0(W \nW) itera- 
tions of the PCGM. 

We now return to the steps involved in solving the system of equations (6.2). As a first 
step it would be necessary to solve the much smaller system of equations (6.7). Here the 
dimension of the vector Z B is p, the number of vertices of the domain £2. Now to be able 
to solve (6.7) to an accuracy of 0(e~ bw ) using PCGM the residual 

R B = SU B -h B 

needs to be computed with the same accuracy and in an efficient manner. The bottleneck 
in computing R B consists in computing (Aji)~ 1 AibUb to an accuracy of 0(e~ hw ) and it 
has already been seen that this can be done using 0(W InW) iterations of the PCGM for 
computing (Au)~ 1 AibUb for a given vector U B . 

We now show that it is possible to explicitly construct the Schur Complement matrix § 
in 0(W \nW) iterations of the PCGM. 8 is a p x p matrix. Let e k be a column vector of 
dimension p with 1 in the kth place and elsewhere. Let & k = Se k . 

Then the Schur Complement matrix S can be written as 

S=[Si,S2,...,S p ]. 
Now by a well known result on the Schur Complement we have 
UgSU B = min V T AV 



v. v B =u B 



= min 

h] 1,1 



Here U B = [gi,g2, ■ ■ ■ ,gp] T ■ Hence using Theorem 3.1 we conclude that 



And so we obtain 

IIS" 1 !! <C(lnW) 2 . (6.14) 
Here the norm denoted is the matrix norm induced by the Euclidean norm. Now 

§k = Se k = {A BB -A T IB A u l A IB ) e k . 

Let (§/t) a be the approximation to S k computed using (9(WlnW) iterations of the PCGM 
to compute A^Ai B e k . Then 

||§,-§£|| = 0(e- w ). 

Let § a denote the matrix 



S" = [S?,^,...,S fl p ]. 
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Clearly 

||S-$ fl || =0(e- bw ). 

Now to compute S a requires 0(W In W) iterations of the PCGM since p is a fixed constant. 
Hence we can solve (6.7) as 

§Z B = h B 

by replacing § by the matrix §" . Let Zg be the solution of 

§"Z£ = h B . 

Since 

S" = S + c5S, 

we have 

(ST 1 = (I + S- l 8S)S- 1 . 

Thus 

IIS" 1 - (SVll < 2||§- 1 || 2 ||5S|| < 0((lnW) 4 )||5§|| 

for ||5S|| small enough. 
Hence 

||S -1 -(S fl ) -1 || = 0(e- bw ). 

Therefore 

\\? B -Z B \\=0(e- bw ). 
Having solved for Z B we obtain Zj by solving 
A//Z/ = h[-A IB Z B 

using (9(WlnW) iterations of the PCGM. Hence the solution Z can be obtained to expo- 
nential accuracy using (9(WlnW) iterations of the PCGM. 

We shall now briefly examine the complexity of the solution procedure for the h-p 
finite element method. Since finite elements have to be continuous along the sides of the 
elements, the cardinality of the set of common boundary value is large in the h-p finite 
element method. Let § denote the Schur Complement matrix for the h-p finite element 
method. In H1I10I it has been shown that an approximation S a to S can be obtained such 
that the condition number % of the preconditioned system satisfies 

*<C(l + (lnW) 2 ), 

where C denotes a constant. Then to solve SU B = h B to an accuracy 0(e~ hw ) will require 
(9(WlnW) iterations of the PCGM using §" as a preconditioner. Now to compute the 
residual in the Schur Complement system to an accuracy of 0(e~ bw ) requires 0(W) iter- 
ations of the PCGM to compute AJ^AisVe. Hence we would need to perform 0(W 2 InW) 
iterations of the PCGM for computing AJ^Vi, where V/ will vary after every sequence of 
(9(WlnW) steps. So the h-p finite element method requires (9(W 2 lnW) iterations of the 
PCGM to obtain the solution. 

Hence the proposed method is faster than h-p finite element method by a factor of 
0(W). 
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7. Technical results 

Lemma 1A. Let {{"^(v^,^)},^,^^ 1 (|,tj)}/} e n M > w . Then there exists {{v k tj 
(v k ,fa)}ij.k, K l (£,T?)M such that v\ x = Ofor all i,k, e H 2 (Cl k Lj )for 2 < j <M 
and all i and k, vf +1 EH 2 (S)forl = 1,2,...,L and w = u + veH 2 B ' 2 (Q.). Moreover the 



estimate 



P 4 M L 
III(PMf +1 -0- 2A 1vt(v fe ,^||^ +£||vr 1 (« I i?)fc 

k=\i=\j=2 /=1 

\*= 1 y s c£i*uflj=, J u(r s )<°° 

x(ii[« ft ]ii 2 fe +ii[(<)iii^+ii[K/]ii? /2 , fe ) 



+ E (ii[»^ 1 ]t + ii[« 1 ) fl ]ii? A , + ii[« 1 ) fl ]ii? Alt : 



/ p h M 

\fc=l!=l .7 = 2 hj 

+iH + \^n)\\ 2 2 ] (7.D 

i=i ' / 

/zoWs. i/ere e w is exponentially small in W. 

We first make a correction {{^-(v*, <pk)}i.j,k, {r P t +i (£, ,!?)}/} such that rf j = for all i 
and k and at the vertices Qi for / = 1, ... ,4 of £2y, 

+ =*(&). 
(("yK + (^)vJ(a) = "v,(a), 

+ (4')fc)(e/) = (7.2a) 

provided g/ is not a vertex of for all i,k. If 2/ is a vertex of choose rf 2 such that 

(«" 2 +4)(e/) = »n(G/), 
((^2)v,+(^ 2 )vj(a) = (4ik(a), 

(("* 2 )ft + - («n)fc(e/)- (7.2b) 

Here s(<2/) denotes the average of the values of s at 2/ over all the elements which have 
Qi as a vertex. 

We can find a polynomial rfAVk^k) on such that r k -(Qi) = ai,(rfj) Vk (Qi) = 
bi,(r k j)fa(Qi) — ci for I = I,. . . ,4. Here the values a/,fc/,c; are defined by (7.2). More- 
over is a polynomial of degree less than or equal to four and the estimate 
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K,-(v*,fe)||* <c(f>/| 2 + N 2 +k,| 2 ) (7.3) 



2,n ; . ■ i , , 
hi \l=i 



holds for j > 2 and all i and fe. Next consider £lf +1 S fi p+1 . Now 

K +1 ), 2 = («f +1 )^, 2 + («f +1 ) I) 7 7;t2 . 

Let ^ p ^jt 2 , r) Xl and f) T2 denote the polynomials of degree W in | and 77 separately which 
are the approximations to ^ p £> Xl , r\ X[ and T7 r , in the space of polynomial of degree W as 
defined in Theorem 4.46 of |[T2l . 

Let Pj,j=l,...,4 denote the vertices of S. Then | X( (P ; ) = E, Xj (Pj) and r) x , (P/) = T] Xi (P/) 
for i = 1,2 and 7 = 1 , . . . , 4. Now 

( M f +l )« =(uf +1 )^ xi +(uf + \n xi ,^d 

Hence (uf + )".(P/) = {uf + )x, -(Pj)i f° r ' = 1,2 and jf = 1, ... ,4. Therefore we can find a 
polynomial r( +1 (%, 77) onS= (Mf +1 )" 1 (i2f +1 ) such thatfor y = 1,...,4, 

K +1 +rf +1 )(P 7 )=«( J P,), 

(( M f +1 ) JCl + (rf +1 ) JCl )(P ; -)=«, 1 (/' ; -),and 

((«f +1 ) X2 + (rf +1 ) X2 )(P,)-«, 2 (^)- 
Now let q(y) be a polynomial of degree W defined on / = (—1,1). Then by Theorem 4.79 

of ma 

\\ q \\l {I) <C(\ n W)\\q\\l 2i . (7.4) 
Here C is a constant. Hence using (7.3) and (7.4) we obtain 

p h M L 

III(PMf +1 - ; )-^ll^(v fc ,^)||^ +IIK +1 (§,T7)|| 2 2 

k=\i=l j=2 ' ' ' 1=1 

<K(lnW)(£ £ rf(A,,y s )- 2 ^ 
V fc=1 y s c£i*UB^, J u(y s )<o<. 

xai^i^+iiK^Jii^+iiK^m^) 

+ E (iiK +1 ]ii 2 yi + ii[(< +1 )iii 2 /2 . 7l + ii[K +1 r]ii 2 /2ri )) 
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( p l k M 

+ £ W III(PMf +1 ^)- 2 ^ll"t-(v*,«-^|P 

\*=ii=l;=2 



Let 



(7.5) 



i=\ 



ytj(Vk, <h) = u tj(Vk, fa) + rfjiVk, fa) and 
yf +l (^,i 1 )= U f +l ^ 1 i 1 ) + rf +l (^ 1 ). 
Now we define a correction {{s^ ;(Vk,fa)}ij t ki {*/ + (^i 1 ?)}/} sucn that s n = for all 



i and sfj G H 2 (£2fj) for 2 < ; < M and all z and /t, sf" 1 e H 2 (S) for / = 1 , . . . ,L and 
w = y + s£H 2 fj 2 (£1). 

Consider £2* , with 2 < j <M. Let 

1 



Hfa) 



Mj-yfj + i)\ 



gm) = -jWj-yij+OvJf, . and 



Hi{<l>k) = -^(yij-yij+ih k \ r 



(7.6) 



In the same way we define F/ , G/ and /// for I = 1, . . . , 4. If Ji C <9£2 for some l,Fi,Gi and 
/// are defined to be identically zero on //. Now F/, G/ and Hi are polynomials of degree 
W that vanish at the end points Qi and Qi + 1 of fi ■ If 73 C <9i2^ j n d£l k i2 for some /. A: then 
the factor of | will be missing in the definition of F^(fa), Gi(fa) and H^fa). We wish 

ft for 



to define sfj(v*,&) on &y such that = l 



: G, and 



/ = 1,...,4. 
We 

A=-D|"| 5=0 , g/=D(j«| JJ=0 fori* ED{ 



We now cite Theorem 1.5.2.4 of 0. The mapping u — > {{/it}™ = o > {£*}*L=o } defined by 

< K + ) has a unique continuous extension as 



A 



■t 



3 

A 

Q 2 



v 

k 



Q =(«,P) 
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an operator from W]"(R + x R + ) onto the subspace of 

m— 1 m— 1 

F = f[ W m -^ 1/p (M + ) x ]~[ W m -'- 1/p (K + ) 
k=0 1=0 

defined by 

(a) D l n f k (0) = D^gi (0) , I + k < m - 2/p for all p ^ 2, and 

(b) J 5 - D|g,(f ) | 2 df /f < oo, / + A; = m - 1 for p = 2. 

Hence using a partition of unity argument it is enough to show that 

(i) J(f \D h F x {t + P) -H 2 (a-f)| 2 dfA, and 

(ii) J S |Gi(r + J3) -£> Vjt F 2 (a-f)| 2 d*A. are nnite - 



Conditions (i) and (ii) follow by applying the above theorem to a neighbourhood of the 
vertex 2j = (a,jS) of£2f ; .. 
Now 

j( |D fcJ F 1 (f + j3)-// 2 (a-f)| 2 df/f 

<2 f 8 \D h Fi(t + p)\ 2 dt/t + 2 [ S \H 2 {a-t)\ 2 dt/t. 



Moreover from Theorem 4.82 in 01211 we have that if q(y) is a polynomial of degree IV 
on/= (—1,1) such thatg(-l) =q(l) — 0, then 



2 



Now by (7.4), 



Ml {1) <KlnW\\ q \\l 1/2(i . 
Hence we conclude that 

« [ p^t + P) -H 2 (a-t)\ 2 dt/t < C(lnW) 2 (||D fe F 1 ||2 /2ft + im] 2 ^. 



(7.7) 



A similar result holds for (ii). 

Hence we can define {{sf /(v*, (j>k)}ij.k, (4 , 1 7)}/} such that *f i = for all i and £, 
si j e H 2 (D.lj) for ; > 2, sf +1 e // 2 (5) and w = y + s E fl^(Q). 

Let v^V*, fc) = rfj(v k , fc) + s*/^, 0,) and vf +1 , 77) = rf +1 (| , 17) + *f +1 ,Tj). 
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Now from (7.7) we conclude that there is a constant K such that 

P h M 



p+i 



(4,1)11 



t=n=i j=2 



Ui l=l 



i 

2,S 



<K(\nWf £ £ d(A h% )- Uk 

\fe=ly s cn*u4,M(y J )<=o 

xdi^lilL + |i[(«* w " 2 1 w " 2 



*V J 11 1/2,* 



+ I (ii[« p+1 ]iiL + ii[K +1 r]ii? /2 , R + ii[(< +1 niiw,J 



y s C£2' 



*2 I >"l/2,y s > 



( P 4 M 

+ £ w III(pMf +w r n *iK.-(^«-% 112 



\yt= 1 1= 1 j=2 



-I K +1 (4, dii 
/=i 



2 

2,5 



(7.8) 



Combining (7.5) and (7.8) gives the estimate (7.1). 
We now prove the last result of this section. 

Lemma 7.2. Let w = u + v e H 2 f (Q) . Here {{«* AVki§k)}i {u'', +> (4,1])}/} £ n"'"' and 
{{v^(V£, (j>k)}ij.ki { v 1 + (4 1 1 ?)};} & ai defined in Lemma 7.1. 77zen f/ze estimate 

2 



Ml 2 3, 



ff 2 ' 2 (r[°l) 



<9w 

aw 



Hf2 (rW) 



<C(lnW) 2 ( £ N 2 +E £ £ 

\*: an*nr[o]^a /€^t=/-l r5 can*nr /:1 u(i> s )<~ 

x d (A t , %) - 2A * ( 1 1 «* - h k 1 1 § A + 1 1 < 1 1 \, 2A ) 



II I ^(A fe ,r,)- 2 ^ 

l€jYk=l-\ y s CdQ. k r\T,,p.(y s )<<x> 



+ 1 I 



P+l||2 



+ 1 I 

/e ^ rs c5n p+1 nr ; 



fdu" 



+1 



dn i A i 



1/2*/ 



1/2,* 



+ 1 I 
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x d(A k , Ys y 2 \\\[u k ]\\l ?s + IIKOlll^ + \\H k n\ 2 lM ) 

+ E (iiK +1 ]iio^ + ii[(< +1 nii^ + ii[K +1 niiw,J 



(p p h M 

E N 2 + E E E(PMf +1 ^r 2A i«f, 7 (v„^) -m 2 
t=l k=\i=\j=2 'hi 

+ thl + \$*)\\l)j (7-9) 

/zoZcfc. Here e w is exponentially small in W. Now 
IMI 33 = inf {\\q\\ 2 ,2 (0 .}■ 

„ 2 2 (r [0] } «l r[0 ]=W ^ (O) 

Let G C 2 (K) such that t = 1 for r k < pfi k and k = for r<. > p. Let g/, = gfy and 
^0 = 1- LLi Let 9o = 1 - L[' = i Define = {x: d{A k ,x) < pn k } fovk=l,...,p 
and let Q. p+l = Q. \ (J k=l & k • Then it can be concluded that 



IN 2 



sf'2(rW) 



<c[ E inf {ll*H 2 22 k }+ inf {IM 2 } 



(7.10) 



Now using (3.9) we have 



IN| 2 2 2 t ^Cd^p+llfe^^-^e-^-ft^H 2 ). (7.11) 

P 

Let us choose the cut-off function 9 k to be a piecewise polynomial such that 
9k(Vk) = 1 for v k < ln(p^), 
%(ln(pM*)) = l,^ (1) (ln(p/Ht)) =O,0f(ln(p M *)) =0, 
fc (lnp) = 0, fc (1) (lnp) = 0, 0f '(lnp) = 0, and 
0k(v k )=0 for v k >lnp. 

Here 0^ denotes the Zth derivative of k with respect to v k . Then 0^ is a polynomial of 
degree five in v k forln(p/4) < v k < lnp. Now using (7.10) and (7.11) we have 
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E inf {||ft|| 2 „ , } 



=,v% '"V' 2 (a*) J 



<C\\ £ \h k \ 2 + £ \\(q k (v k M)-h k )e- 2 ^% u _ 



k: r^nrPl/B 



+ £ llfe(v*,^)-^)e- 2(1 - &)Vt ll 



3/2,(-~,tap) 



/t: r t+1 nrl°]^0 

Let 77* = In p + (M + 1 - ;) In jx k and l) = {r\)_ x , rf). Then 



\\(q k (v k y u )-h k ) e - 2 ^-^\\ 



M+l 



<C E ll(ft(v*,^)-A t )e 

I .7=2 



3/2,(-~,fap) 



M+l 



+ EIK(^ ; ^)-Me- 2(1 ^K" 2 



.7=2 



+ E/ ^(fe(^^)-^)e^ ft) ")(r7} + cT) 

7—2 



((^(,, V 4)-^) e - (1 - ft)S )(T ] }- ( 7) 



ds 



da 



Here 5 > 0. Now 



Jo 



ds 



((q k (s,ri)-h k )c-^)(rf + cj) 



-^((^(^^)-^)e- (1 - ft)s )(T?*-<T) 



da 
a 



<^(E/||rfe-%)(^ + -^) 



d' 



^(9*-A*)fa*-ff>Ytf) 



' e -2(i-ft)(^+d)da 
a 



1 J' 

/=0 J0 



x ( e -2(l-&)(ijJ+a) _ e -2(l-A)(ij)-<r)jdCT^ 
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Hence 

\\(q k (v k M)-h k )e- 2 ^Xn,-. 



3/2,(— MP) 

—IXi, 



>y s crP]nr t+1)/i (f s )<oo 



x[\\(q k -h k )\\l fs - 



dv k 



2 

Till vv || \qk-' l k)\\ 2 

1/2,15 



I . (7-13) 



This follows from Theorem 4.82 in 1 12 1 which states that if p(y) is a polynomial of degree 
Niny such that p(l) = p{— 1) = 0, then 

- 1 P 2 (3'),../^ 1 „, rll „„2 



/ f^4dy<Cln^V||p|| 
J -i 1 - y z 



-1 1 -y 2 J ~ "i-W 

Now ^(Vjt, y£) - % = fyw^, y£) -% = ^(v*, y/£) - %) + (0* - 1)% 

Moreover w = r + s + u, as has been defined in Lemma 7.1. Here w(v£, i//^) is a poly- 
nomial of degree W, s(v^, l//*) = and r(Vfc, ) is a polynomial degree four. 
Hence using (7.4) and (7.13) we conclude that 



<c 



£ d(A ft , 7 ,)-^((lnW) 2 (||( M (v„^)-^)||^ A 



+ |/ ! ,| 2 )+lnW||r(v,,^)||2 /2 .)). (7.14) 
Hence using (7.5) and (7.14) it can be concluded that 
E inf {||ft|| 2 2 2 t } 



< C(lnW) 2 ( £ \h k \ 2 + E E E «*(4t> 

x(ll"*-^ll 2 A + ll<ll? /2A ) + E E d^y,)- 2 ** 

£=iftcn i -uB*,ji(j',)<°o 

x(ii[^]ii^ + ii[(<jiii 2 /2 .+ii[(<r]ii 2 /2 .) 



E (ll^ll^ + ll^nil^ + IIK^ 1 ) ]!!^) 
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P M h 



yk=l 7=2 1= 1 



2,Q K . 



+en 2 +ek +1 (^)iil 



1=1 



(7.15) 



In the same way we can show that 



Now 



Here 



inf 



<C(lnW) 2 £ N 2 +£ £ £ 

\e a£i*nrlo]/B ie®k=i-i rs cdn k nr h n(%)<<*> 

x^^-^dl^-^ilo^+lKiiWJ+I E 



xJ(A,,7 s )- n HII[" ft ]ll 2 fe + ll[«r]ll? /2 , fe + 



+ E (iiK +1 ]ii ^ + ii[(< +1 )iii 2 /2 , r! + ii[K +1 r]ii 2 /2 , r! ) 



y s C£l 



E E 



p+l||2 



' e ® ri can" +1 nr ; 
p m 4 

EEE 

yt=l 7=2 i=l 
L 

E 

*:=i i=i 



duP+ 



1/2,7*/ 



+<v ( E E E(p< +w r 2A *ii<An,<y -m 2 



+ EN 2 +EK +1 (^)l 



2 

2,5 



(7.16) 



<9w 
dN 



inf 

»l r [i] = (^), 



112 +Lii*/»^< m - 



"l2 (£J) 



|o|=l 
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Let 9k £ C 2 (R) be as defined earlier and q k = q 9 k . Let qo = 1 — Then, as before 



«J' 3 (ri 1 !) 



<C 



inf 



) + 



inf 



* l a£2*nr[ 1 ] = (^)/* "P ^ <?ol an/'+inr[i] = (^ 



Now 



(£1*) 



Jo* «„ V 

for J3jt > 0. And so 

[ e 2 P« v *\q k (v k ,<k)\ 2 dv k d<l> k < K \\q k \\ 2 



(£1*) 



Hence for < fi k < 1, there exists a constant C such that 



inf 



<E 



inf , \ (ll*H 2 n k ) 



< C £ inf 



( e2ftV * 1 1 [ P$U«*( V *' <M | 2 dV fc d^ ^ 



Thus by similar arguments as before it can be shown that 

2 



fdw\ 
\dNj A 



<C(lnW) 2 ( EE E d(At, 
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fdu k " 


a 




I 



E E 

d(A k ,y s )- n * 



dN 



+ E E 



x(iiK]ii^ + ii[(<r]ii^ + ii[(<niik) 



l/2,7s 



+ E (ii[" p+1 ]iiL + iiK +1 r]iii 2ri + i![K +1 r]iii 2ri ) 

y s coF +l 



( P M h 

V= l./=2 i=l ' ,J 



+ EN 2 +EK +1 (<^)ii2 



p 

E 

>t= i /=i 

Combining (7.15)— (7.17) we obtain the required result. 



(7.17) 
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